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We describe some connections, via composition, between two functional spaces: the 
space of (sub)critical branching mechanisms and the space of Bernstein functions. The 
functions : a; — > j;" where x > and < a < 1/2, and in particular the critical 
parameter a = 1/2, play a distinguished role. 

1 Introduction 

This note is a prolongation of [8] where the following remarkable property of the function 
ec : X — > a;" was pointed at for a = 1/2: if \E' is a (sub)critical branching mechanism, then 
\l/ o ei/2 is a Bernstein function (see the next section for the definition of these notions). In 
the present work, we first show that this property extends to every a g] 0,1/2]. Then we 
characterize the class of so-called internal functions, i.e. that of Bernstein functions $ such that 
the compound function o $ is again a Bernstein function for every (sub) critical branching 
mechanism ^. In the final section, we gather classical results on transformations of completely 
monotone functions, Bernstein functions and (sub) critical branching mechanisms which are 
used in our analysis. 

2 Some functional spaces 

2.1 Completely monotone functions 

For every Radon measure on [0, oo[, we associate the function :]0, oo[^ [0, oo] defined by 
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i.e. is the Laplace transform of /i. We denote by 



CM := {C^ : C^{q) < oo for all g > 0} , 



(2) 
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which is an algebraic convex cone (i.e. a convex cone which is further stable under inner 
product). The celebrated theorem of Bernstein (see for instance Theorem 3.8.13 in [6]) identifies 
CM with the space of completely monotone functions, i.e. functions / :]0, oo[— > [0, oo[ of class 
such that for every integer n > 1, the n-th derivative /^"-' of / has the same sign as (—1)". 
Recall from monotone convergence that has a (possibly infinite) limit at 0+ which coincides 
with the total mass of //. 

We shall focus on two natural sub-cones of CM: 

Bi:^\c.:[ (lAx-^)ii(dx)<oo\ (3) 
[ J]oM J 

We further denote by B| the sub-space of functions in Bi which are the Laplace transforms of 
absolutely continuous measures with a decreasing density : 

B]; := : fj,{dx) = g{x)dx,g decreasing and J (1 A x~^)g{x)dx < oo| . (4) 

Note that the density g then has limit at infinity. 



2.2 Bernstein functions 

For every triple (a, b, A) with a,b> and A a positive measure on ]0, oo[ such that 

/ (x A l)A{dx) < oo , (5) 

J]Q,oo[ 

we associate the function $a,6,A :]0, oo[— > [0, oo[ defined by 

*aAA(5) ■.= a + bq+ f (1 - e-'i-)A{dx) , (6) 

J]0,oo[ 

and call $a,6.A the Bernstein function with characteristics (a, b, A). We denote the convex cone 
of Bernstein functions by 

B2 := {^a,b,A ■ a,b>0 and A positive measure fulfiUing (5)} . (7) 

It is well-known that B2 can be identified with the space of real-valued C°° functions / : 
]0,oo[— > [0,oo[ such that for every integer n > 1, the n-th derivative /*^"^ of / has the same 
sign as (—1)"^^. Sec Definition 3.9.1 and Theorem 3.9.4 in [6]. 

Bernstein functions appear as Laplace exponents of subordinators, see e.g. Chapter 1 in 
[3], Chapter 6 in [9], or Section 3.9 in [6]. This means that $ G B2 if and only if there exists 
an increasing process a = {(Jt,t > 0) with values in [0, 00] (00 serves as absorbing state) with 
independent and stationary increments as long as at < 00, such that for every t > 

E(cxp(-g(T()) = exp(-t$(g)) , g > 0. 

In this setting, a is known as the killing rate, b as the drift coefficient, and A as the Levy 
measure. 

We shall further denote by B2 the subspace of Bernstein functions for which the Levy measure 
is absolutely continuous with a monotone decreasing density, viz. 

B2 := |$a,fe,A : a, 6 > and A{dx) = g{x)dx, g >0 decreasing and J {x A l)g{x)dx < ooj . 
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2.3 (Sub)critical branching mechanisms 



For every triple (a, b, U) with a,b> and U positive measure on ]0, oo[ such that 

/ {xA x^)U{dx) < oo (8) 

^]0,oo[ 

we associate the function ^a,6,n :]0, oo[^ [0, oo[ defined by 

*a,6,n(g) aq + bq''+ [ (e"^^ - 1 + qx)U{dx) , (9) 

and denote the convex cone of such functions by 

B3 := {^a,b,n ■ a,b >0 and 11 a positive measure such that (8) holds} (10) 

Functions in B3 are convex increasing functions of class C°° that vanish at 0; they coincide 
with the class of branching mechanisms for (sub)critical continuous state branching processes, 
where (sub)critical means critical or sub-critical. See Le Gall [7] on page 132. 

Alternatively, functions in the space B3 can also be viewed as Laplace exponents of Levy 
processes with no positive jumps that do not drift to —00 (or, equivalently, with nonnegative 
mean). In this setting, a is the drift coefficient, 2b the Gaussian coefficient, and 11 the image 
of the Levy measure by the map x — > —x. See e.g. Ghapter VII in [2]. 

3 Composition with Gq, 

Stable subordinators correspond to a remarkable one-parameter family of Bernstein functions 
denoted here by (e^^, < a < 1), where 

ea{q) := q"" = " , /"(I - e-^x-'-^^dx , q>0. 
I [L — a) Jo 

Theorem 1 The following assertions are equivalent: 

(i) «G]0,l/2]. 

(n) For every ^ e B3, ^ o G B2. 

The implication (ii) =^ (i) is immediate. Indeed, ^0,1,0 '■ Q ^ belongs to B3, but — 
^0,1,0 is in B2 if and only if 2q; < 1. However, the converse (i) ^ (ii) is not straightforward 
and relies on the following technical lemma, which appears as Lemma VI. 1.2 in [8]. Here, for 
the sake of completeness, we provide a proof. 

Lemma 2 For a e]0, 1/2], let o"^"^ = {a^\x > 0) be a stable subordinator with index a with 
Laplace transform 

E (exp (-gcr^"^)) = exp(-a;g") , x,q> Q. 
Denote by p^"^ (x, t) the density of the law of a'^^ . Then for every ,x,t>0, we have 

f) < - r/~(^+'*) 
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Remcirk : The bound in Lemma 2 is sharp, as it is well-known that for any < a < 1 and 
each fixed t > 

i (1 — a) 

More precisely, there is a series representation of p*^")(a;, t), see Formula (2.4.7) on page 90 in 
Zolotarev [10]: 

-na—1 



= - E(-l) r(n + 1)) ^^"(^^")^ 



Using the identity 

r(Q;)r(i - a) 



sm avr 



this agrees of course with the above estimate. It is interesting to note that the second leading 
term in the expansion. 



sin(27rQ;)a; 



27r 

is negative for a < 1/2, but positive for a > 1/2. So the bound in Lemma 2 would fail for 
a > 1/2. 

Proof: In the case a = 1/2, there is an explicit expression for the density 

from which the claim is obvious (recall that r(l/2) = \^)- 
In the case o; < 1/2, we start from the identity 

exp(-xg") = / e-'''p^''\x,t)dt, 
Jo 

and take the derivative in the variable q to get 

roo f 

ag"-^ exp(-a;g") = / e-'^^-p^'^\x,t)dt , 
Jo X 



and then 



a?"-^ (1 - exp(-xg°)) = / e'^* — -r" - -p(°)(x, t)\dt. 

Jo \r(l — a) X ) 



Denote the left hand-side by g{x^ q), and take the derivative in the variable x. We obtain 

roo 

aq^""-^ e-'''p^''\x,t)dt. 
Jo 



dx 

On the other hand, since 1 — 2q; > 0, 



1 roo 
2a-l ^ ^ / qs -2a^ 

^ r(l - 2a) Jo 



and hence 

dx r{l-2a)Jo s^'^Jo ^ y ^ J 
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The change of variables u = t + s yields 

dg{x, q) a r°° ds 



dx r(l - 2a) 

and since g{0,t) = 0, we finally obtain the identity 

^ a „, t 



Jo \r(l — a) X J 

pu PX PCX) rlq 

— dy du ^e-V")(x,ii-s). 

I -2a) Jo ^ Jo io ^ ^ ' ^ 



r 

Inverting the Laplace transform, we conclude that 

a t a r . /■* (is 



r(l-a) a;"^ ' ' ' r(l - 2a) 

which entails our claim. 



We are now able to prove Theorem 1. 

Proof: Let ^a.6,n G B3. Since both aSa and 6e2Q- are Bernstein functions, there is no loss of 
generality in assuming that a = b = 0. Set for t > 

It follows from Lemma 2 that Ua{t) > 0. We have for every q > 

(1 - c-^*)z/«(t)rft 

= / ^{dx)xW--^- 

Jo \ X X J 

= ^o,o,n(ea(g)) • 

As this quantity is finite for every g > 0, this shows that ^o,o,n o e B2. ■ 

Remark : The proof gives a stronger result than that stated in Theorem 1. Indeed, we 
specified the Levy measure z/„ of ^o,o.n ° ^a- Furthermore, in the case a = 1/2, this expression 
shows that ^'o,o,n°6i/2 G B2. It is interesting to combine this observation with the forthcoming 
Proposition 7 : for every ^ e B3, \t' o ei/2 e B2, thus Id x o ei/2) : g ^ Q^iy/Q) is again 
in B3, and in turn ei/2 x o ei/4) e B^. More generally, we have by iteration that for every 
integer n 

62-21-" X o e2-n) e B3 , 

and 

ei_2-n X o e2-n-i) e B2 . 
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4 Internal functions 



It is well-known that the cone CM of completely monotone functions and the cone B2 of 
Bernstein functions are both stable by right composition with a Bernstein function; see Propo- 
sition 8 below. Theorem 1 incites us to consider also compositions of (sub)critical branching 
mechanisms and Bernstein functions; we make the following definition : 

Definition 3 A Bernstein function $ e B2 is said internal i/^^ o $ e B2 for every e B3. 

Theorem 1 shows that the functions are internal if and only if a g]0, 1/2]. The critical 
parameter a = 1/2 plays a distinguished role. Indeed, we could also prove Theorem 1 using the 
following alternative route. First, we check that ei/2 is internal (see [8]), and then we deduce by 
subordination that for every a < 1/2 that o Bq, = o ei/2 o e2a is again a Bernstein function 
for every e B3. Developing this argument, we easily arrive at the following characterization 
of internal functions : 

Theorem 4 Let $ = $a,6,A € B2 6e a Bernstein function. The following assertions are then 

equivalent: 

(i) $ is internal, 

(ii) $2 G B2, 

(iii) 6 = and there exists a suhordinator a = {at,t > 0) such that 

roo 

A(dx) = c / t'^^^Fiat e dx)dt . 
Jo 

Proof: (i) ^ (ii) is obvious as ^0,1,0 o $ = 

(ii) =^ (i). We know from Theorem 1 or [8] that for every ^ e B3, ^ o ei/2 G B2. It follows 
by subordination that for every Bernstein function k, G B2, ^ o ei/2 one B2. Take k, — so 
ei/2 o K — and hence $ is internal. 

(iii) (ii) Let k, denote the Bernstein function of a. We have 

$(g) = a+ f (1 - e-'^'^)A{dx) 

J]0,oo[ 

^ a + c [ r dt(l - e-nt-^^^F(at G dx) 

J]0,QO[JO 

= a + c rdi(l-e-*''(^))r3/^ 
^0 

The change of variables tK{q) = u yields 

$(g) = a + c'y^<^ 

and hence 

^^(q) = + 2ac'^/K{q) + c''^K{q) . 
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Since t^l"^ = ei/2 o is again a Bernstein function, we thus see that e B2. 
(ii) ^ (iii) Recall that the drift coefficient h of <l>a,6,A is given by 

lim $a,b,A(g)/g = 6; 

see e.g. page 7 in [3]. It follows immediately that 6 = whenever := $a5A ^ -^2- Recall from 
Sato [9] on page 197-8 that if r^^) and r*-^-* are two independent subordinators with respective 
Bernstein functions ^'^^^ and ^'^^^ then the compound process t^"^^ o r^^'' := r*^^-* is again a 
subordinator with Bernstein function := o moreover its Levy measure 

A(3) 

IS 

given by 

A^^\dx) = /"Pfr/'^ e dx)A^^\dt) , 
Jo 

where A^^^ denotes the Levy measure of t^'^\ As $a,6,A = ©1/2 o k, and the Levy measure of ei/2 
is ct'^/'^dt with c = 1/(2-^/7?), we deduce that 

/•oo 

Aidx) = c / P(cTt e dx)r^'^dt . 
Jo 

The proof of Theorem 4 is now complete. ■ 



It is noteworthy that if $a,b,A is internal and A ^ 0, then 

xA{dx) = 00 . 



/ 

^10. 



00 . 



/]0,oo 

Indeed, 

/ xA{dx) ^ c r f xF{at e dx)r^''^dt = c f°° E{ai)r^^'^dt = 

J]0,oo Jo J]0,oo[ Jo 

For instance, the Bernstein function q log(l + q) of the gamma subordinator is not internal. 

Corollary 5 For every G B3, wewrite $ for the inverse function of ^ and then for its 
derivative. Then 1/$' is internal. 

Proof: It is known (see Corollary 10 below) that 1/$' is a Bernstein function; let us check 

that its square is also a Bernstein function. 

We know that ^" G Bi (Proposition 6 below) and $ G B2 (Proposition 9 below); we 
deduce from Proposition 8 that o $ g Bi. If we write /(/) : x — > Jq f{y)dy for every 
locally integrable function /, then again by Proposition 6, we get that o $) is a Bernstein 
function. 



Now 



so 



and we conclude that 



o $ = 



$' O ^f)3 ' 



($')3 

1 



2($')^ 



/(^" o $) G B2 . 
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5 Some classical results and their consequences 



For convenience, this section gathers some classical transformations involving B^, j e {1,2,3} 
and related subspaces, which have been used in the preceding section. We start by considering 
derivatives and indefinite integrals. The following statement is immediate. 

Proposition 6 Let j = 2,3 and f :]0, oo[— > [0, oo[ be a C°°-function with derivative f. For 
j — 3, we further suppose that limg_>o f{q) — 0. There is the equivalence 

feBj ^ /'eB,_i. 

The next statement is easily checked using integration by parts. 

Proposition 7 Let j = 2,3 and consider two functions f,g :]0, oo[— >• [0, cx3[ which are related 
by the identity f{q) = qg{q). Then there is the equivalence 

f e Bj and lim /(g) = <^ y e B]_i . 

Proposition 7 has well-known probabilistic interpretations. First, let cr be a subordinator with 
Bernstein function / € B2 with unit mean, viz. E(o"i) = 1, which is equivalent to /'(0+) = 1. 
Then the completely monotone function g{q) := f{q)/q is the Laplace transform of a probability 
measure on M+. The latter appears in the renewal theorem for subordinators (see e.g. [4]); 
in particular it describes the weak hmit of the so-called age process A{t) — t — gt as t ^ 00, 
where gt :— sup {cr^ : as < t} . Second, let X be a Levy process with no positive jumps and 
Laplace exponent / G B3. The Levy process reflected at its infimum, — info<s<f X^, is 
Markovian; and if r denotes its inverse local time at 0, then cr = —X o r is a subordinator 
called the descending ladder-height process. The Bernstein function of the latter is then given 
by g{q) = f{q)/q'i ; see e.g. Theorem Vn.4(ii) in [2]. 

We next turn our attention to composition of functions; here are some classical properties 

Proposition 8 Consider two functions f, g :]0, oo[— > [0, oo[. Then we have the implications 

f,g eB2 ^ fogeB2, 

f e CM andgeB2 =^ fogE CM , 
f eBi andgeB2 =^ / o ^ e Bi . 

The first statement in Proposition 8 is related to the celebrated subordination of Bochner (see, 
e.g. Section 3.9 in [6] or Chapter 6 in [9]); more precisely if a and r are two independent 
subordinators with respective Bernstein functions f^ and f-r, then uor is again a subordinator 
whose Bernstein function is f^. o f^. The second statement is a classical result which can be 
found as Criterion 2 on page 441 in Feller [5]; it is also related to Bochner's subordination. 

Finally we turn our attention to inverses. 
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Proposition 9 Consider a function f :]0, oo[— >-]0, oo[. Then 

/eB2UB3 =^ l//eCM. 

Further, if denotes the inverse of f when the latter is a hijection, then 

feBsJ^O =^ f-'eB2. 

We mention that if / G B3, the completely monotone function 1// is the Laplace transform 
of the so-called scale function of the Levy process X with no positive jumps which has Laplace 
exponent /. See Theorem VILS in [2]. On the other hand, f~^ is the Bernstein function of the 
subordinator of first-passage times Tt :— inf {s > : > t}; see e.g. Theorem VILl in [2]. 
Finally, in the case when / G B2 is a Bernstein function, the completely monotone function 
1// is the Laplace transform of the renewal measure U{dx) — f^f[(Tt G dx)dt, where cr is a 
subordinator with Bernstein function /. 

CorollEiry 10 Let ^ ^ be a function in B3, and denote by ^ — G B2 its inverse 
bijection. Then q — > l/^'{q) and Id/^ : q — > q/^{q) are Bernstein functions. Furthermore 
: q — > l/($(g)$'(g)) is completely monotone. 

Proof: We know from Propositions 6 and 9 that both $ and are Bernstein functions. We 
conclude from Proposition 8 that 1/$' = ^' o $ is again in B2. 

Similarly, we know from Proposition 7 that g — > is a Bernstein function, and com- 

position on the right by the Bernstein function $ yields Id/$ that is again in B2. 

Finally, we can write !/($$') = / o $ where f{q) = ^'{q)/q. We know from Proposition 6 
that G B2, so / G CM by Proposition 7. Since $ G B2, we conclude from Proposition 8 that 
/o$gCM. ■ 



If $ = is the Bernstein function given by the inverse of a function G B3, the Bernstein 
function 1/$' is the exponent of the subordinator defined as the inverse of the local time 
at of the Levy process with no positive jumps and Laplace exponent \E'. See e.g. Exercise 
VIL2 in [2]. On the other hand, Id/$ is then the Bernstein function of the decreasing ladder 
times, see Theorem Vn.4(ii) in [2]. The interested reader is also referred to [1] for further 
factorizations for Bernstein functions which arise naturally for Levy processes with no positive 
jumps, and their probabilistic interpretations. 

Next, recall that a function / :]0, oo[— > R+ is called a Stieltjes transform if it can be expressed 
in the form 

/(,) = 6+/ ,>0, 

i[o,oo[ t + q 

where 6 > and u is a, Radon measure on R+ such that /[o,oo[(l ^ t~^)iy{dt) < 00. Equivalently, 
a Stieltjes transform is the Laplace transform of a Radon measure /i on of the type ii{dx) — 
bSo{dx) + h{x)dx, where 6 > and h is a completely monotone function which belongs to 
L^[e~'^^dx) for every q > 0; see e.g. Section 3.8 in [6]. 
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Corollciry 11 Let / e B2 6e a Bernstein function such that its derivative f is a Stieltjes 
transform. Then for every Bernstein function g e B2, the function fo^ is completely monotone. 

Proof: We can write 

rq roo 

f(q) ^ a + bq+ dr dxe~'''^h(x) , ? > , 
Jo Jo 

where a, 6 > and /i e Bi. Thus 

f(q) = a + bq+ r dx(l-e-'^'')^, q>0, 
Jo X 

and then 

/ o -{q) r - e-^/^(«))^ . 

9 9{q) Jo X 

We already know from Proposition 9 that a + b/g G CM. The change of variable y = x/g{q) 
yields 

y 

For each fixed y > yg is a. Bernstein function, so by Proposition 8, the function q — > h{yg{q)) 
is completely monotone. 

We conclude that for every integer n > 0, 
which establishes our claim. ■ 



r dx{l-e-^'3^^^)^ = r {I - e~y)h{yg{q))- 
Jo X Jo 
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